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Abstract. The tropical variety of a d-dimensional prime ideal in a 
polynomial ring with complex coefficients is a pure d- dimensional poly- 
hedral fan. This fan is shown to be connected in codimension one. We 
present algorithmic tools for computing the tropical variety, and we dis- 
cuss our implementation of these tools in the Grobner fan software Gf an. 
Every ideal is shown to have a finite tropical basis, and a sharp lower 
bound is given for the size of a tropical basis for an ideal of linear forms. 



1. Introduction 

Every ideal in a polynomial ring with complex coefficients defines a tropi- 
cal variety, which is a polyhedral fan in a real vector space. The objective of 
this paper is to introduce methods for computing this fan, which coincides 
with the "logarithmic limit set" in George Bergman's seminal paper |2j. 

Given any polynomial / E C[x±,X2, ■ ■ ■ ,x n ] and a vector w G M n , the 
initial form in U) (/) is the sum of all terms in / of lowest w- weight; for 
instance, if I = x\ + 22 + %3 + 1 then in( ,o,i)(^) = x i + x i + 1 an d 
in( ,o,-i)(^) = %3- The tropical hypersurface of / is the set 

T(/) = { ro G R n : in„,(/) is not a monomial}. 

Equivalently, T(f) is the union of all codimension one cones in the inner 
normal fan of the Newton polytope of /. Note that T(f) is invariant under 
dilation, so we may specify T(f) by giving its intersection with the unit 
sphere. For the linear polynomial £ above, T{£) is a two-dimensional fan 
with six maximal cones. Its intersection with the 2-sphere is the complete 
graph on the four nodes (1,0,0), (0,1,0), (0,0,1) and -(^, 

A finite intersection of tropical hypersurfaces is a tropical prevariety \12\ . 
If we pick the second linear form £' = x\ + X2 + 2x3 then T{£') is a graph 
with two vertices connected by three edges on the 2-sphere, and T(£)f]T(£') 
consists of three edges of T(£) which are adjacent to — , ^). In 

particular, the tropical prevariety T(£) n T(£') is not a tropical variety. 

Tropical varieties are derived from ideals. Namely, if / is an ideal in 
C[xi, . . . , x n ] then its tropical variety T(I) is the intersection of the tropical 
hypersurfaces T(f) where / runs over all polynomials in I. Theorem 12.91 
below states that every tropical variety is actually a tropical prevariety, i.e., 
the ideal I has a finite generating set /2, . . . , f r } such that 

T(/) = T(/ 1 )nT(/ 2 )n...nT(/ r ). 
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If this holds then {/i, /2, . . . , f r } is called a tropical basis of /. For instance, 
our ideal / = (£, (!) has the tropical basis { x\ +£2 + 2x3, X1+X2 + 2, £3 — 1}, 
and we find that its tropical variety consists of three points on the sphere: 

T(I) = {(1,0,0), (0,1,0), -(-1-1,0)}. 

Our main contribution is a practical algorithm, along with its implemen- 
tation, for computing the tropical variety T{I) from any generating set of 
its ideal I. The emphasis lies on the geometric and algebraic features of this 
computation. We do not address issues of computational complexity, which 
have been studied by Theobald JH]- Our paper is organized as follows. 

In Section 2 we give precise specifications of the algorithmic problems we 
are dealing with, including the computation of a tropical basis. We show 
that a finite tropical basis exists for every ideal /, and we give tight bounds 
on its size for linear ideals, thereby answering the question raised in |16[ 
§5, page 13]. In Section 3 we prove that the tropical variety T(I) of a 
prime ideal I is connected in codimension one. This result is the foundation 
of Algorithm 14.111 for computing T(I). Section 4 also describes methods 
for computing tropical bases and tropical prevarieties. Our algorithms have 
been implemented in the software package Gf an [9, • In Section 5 we compute 
the tropical variety of several non-trivial ideals using Gfan. The tropical 
variety T{I) is a subfan of the Grobner fan of / (defined in Section 2). The 
Grobner fan is generally much more complicated and harder to compute than 
T(I). In Section 6 we compare these two fans, and we exhibit a family of 
curves for which the tropical variety of each member consists of four rays but 
the number of one-dimensional cones in the Grobner fan grows arbitrarily. 

A note on the choice of ground field is in order. In this paper we will 
work with varieties defined over C. In the implementation of our algorithm 
(Section [5J), we have required our polynomials to have rational coefficients, 
but our algorithms do not use any particular properties of Q. It is impor- 
tant, however, that we work over a field of characteristic 0, as our proof of 
correctness uses the Kleiman-Bertini theorem in the proof of Theorem 13. II 

In most papers on tropical algebraic geometry (cf. El El Q22]), 
tropical varieties are defined from polynomials with coefficients in a field 
K with a non-archimedean valuation. These tropical varieties are not fans 
but polyhedral complexes. We close the introduction by illustrating how 
our algorithms can be applied to this situation. Consider the field C(e) of 
rational functions in the unknown e. Then C(e) is a subfield of the alge- 
braically closed field C{{e}} of Puiseux series with real exponents, which is 
an example of a field K as in the above cited papers. Suppose we are given 
an ideal / in C(e)[xi, . . . , x n \. Let /' C C{{e}}[xi, . . . , x n ] be the ideal gen- 
erated by /. The tropical variety T(I'), in the sense of the papers above, 
is a finite polyhedral complex in M n which usually has both bounded and 
unbounded faces. To study this complex, we consider the polynomial ring 
in n + 1 variables, C[e, x\, . . . , x n ] and we let J denote the intersection of / 
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with this subring of C(e)[xf, . . . ,x^]. Generators of J are computed from 
generators of I by clearing denominators and saturating with respect to e. 
The tropical variety of I' is related to the tropical variety of J as follows. 

Lemma 1.1. A vector w G W 1 lies in the polyhedral complex T(I') if and 
only if the vector (l,u>) G M. n+1 lies in the polyhedral fan T(J). 

Thus the tropical variety T(I') equals the restriction of T{J) to the north- 
ern hemisphere of the n-sphere. Note that if / is a prime ideal then so are 
/' and J. Einsiedler, Kapranov and Lind have shown that if I 1 is prime, 
then T(I') is connected. Our connectivity results in Section 3 (which use 
the result of 0) imply the following result which was conjectured in [£]. 

Theorem 1.2. /// is an ideal in C{{e}}[xi, . . . , x n ] whose radical is prime 
of dimension d, then the tropical variety T{I) is a pure d-dimensional poly- 
hedral complex which is connected in codimension one. 

On the algorithmic side, we conclude that the polyhedral complex T(I') 
can be computed by restricting the flip algorithm of Section 4 to maximal 
cones in the fan 7~( J) which intersect the open northern hemisphere in M Tl+1 . 

2. Algorithmic Problems and Tropical Bases 

For all algorithms in this paper we fix the ambient ring to be the polyno- 
mial ring over the complex numbers, C[x] := C[x±, . . . , x n ]. The most basic 
computational problem in tropical geometry is the following: 

Problem 2.1. Given a finite list of polynomials /i, . . . , f r G C[x], compute 
the tropical prevariety T(/i) fl • ■ • fl T(f r ) in W l . 

The geometry of this problem is best understood by considering the New- 
ton polytopes New(/i), . . . , New(/ r ) of the given polynomials. By definition, 
New(/j) is the convex hull in M n of the exponent vectors which appear with 
non-zero coefficient in fa. The tropical hypersurface T{fi) is the (re — 1)- 
skeleton of the inner normal fan of the polytope New(/j). Our problem is to 
intersect these normal fans. The resulting tropical prevariety can be a fairly 
general polyhedral fan. Its maximal cones may have different dimensions. 

The tropical variety of an ideal / in C[x] is the set T(I) := P|j e/ T(/). 
Equivalently, T(I) = {w G M n : in U) (J) does not contain a monomial} where 
in w (J) := (in w (f) : / £ /) is the initial ideal of / with respect to w. Bieri 
and Groves [3] proved that T(I) is a d-dimensional fan when d is the Krull 
dimension of C[x]//. The fan is pure if / is unmixed. In Section 3 we shall 
prove that T(I) is connected in codimension one if / is prime. 

We first note that it suffices to devise algorithms for computing tropical 
varieties of homogeneous ideals. Let h I C Cf xq, x\, . . . , x n ] be the homoge- 
nization of an ideal / in C[x] and h f the homogenization of / G C[x]. 

Lemma 2.2. Fix an ideal I C C[x] and a vector w G W 1 . The initial ideal 
in w (I) contains a monomial if and only if in( 0iU) ) ( I) contains a monomial. 



4 



T. BOGART, A. JENSEN, D. SPEYER, B. STURMFELS, R. THOMAS 



Proof. Suppose x u G in w (I). Then x u = inu,(/) for some / G I. The (0, w)- 
weight of a term in / equals the u>-weight of the corresponding term in /. 
Hence in( 0jtl) )( /) = XqX u G m^ ^( h I) where a is some non-negative integer. 

Conversely, if x u G m^ 0w ^( h I) then x u = in(o iW )(/) for some / G I. 
Substituting xq = 1 in / gives a polynomial in /. The (0, u>)-weight of any 
term in / equals the w- weight of the corresponding term in f\ xo= i. Since 
i n (o,w)(/) is a monomial, only one term in / has minimal (0, u>)-weight. This 
term cannot be canceled during the substitution. Hence it lies in in w (I). □ 

Our main goal in this paper is to solve the following problem. 

Problem 2.3. Given a finite list of homogeneous polynomials f\, . . . , f r G 
C[x], compute the tropical variety T{I) of their ideal / = (/i, . . . , f r ). 

It is important to note that the two problems stated so far are of a funda- 
mentally different nature. Problem l2,ll is a problem of polyhedral geometry. 
It involves only polyhedral computations: no algebraic computations are 
required. Problem 12.31 on the other hand, combines the polyhedral aspect 
with an algebraic one. To solve Problem 12.31 we must perform algebraic 
operations (e.g. Grobner bases) with polynomials. In Problem 12.11 we do 
not assume that the input polynomials /i, . . . , f r are homogeneous as the 
polyhedral computations can be performed easily without this assumption. 

Proposition 2.4. Let I be an ideal in C[x] and let w G W 1 . The following 
are equivalent: 

(1) The ideal I is w -homogeneous; i.e. / is generated by a set S of 
w -homogeneous polynomials, meaning that in w (f) = f for all f G S. 

(2) The initial ideal in w (I) is equal to I. 

Proof. If / has a u>-homogeneous generating set then / C in„,(7). Any 
maximal w-homogeneous component of / G / is in /. In particular in„,(/) G 
/. Conversely, the ideal in,„(/) is generated by w-homogeneous elements by 
definition so, if / = in w (I), then / is generated by w-homogeneous elements. 

□ 

The set of w G R n for which the above equivalent conditions hold is a 
vector subspace of W 1 . Its dimension is called the homogeneity of / and is 
denoted homog(J). This space is contained in every cone of the fan T(l) and 
can be computed from the Newton polytopes of the polynomials that form 
any reduced Grobner basis of /. Passing to the quotient of W 1 modulo that 
subspace and then to a sphere around the origin, T{I) can be represented as 
a polyhedral complex of dimension n — codim(J) — homog(I) — 1 = dim(I) — 
homog(/) — 1. Here codim(J) and dim(J) are the codimension and dimension 
of /. In what follows, T(I) is always presented in this way, and every ideal 
/ is presented by a finite list of generators together with the three numbers 
n, dim(I) and homog(J). 
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Example 2.5. Let I denote the ideal which is generated by the 3 x 3-minors 
of a symmetric 4 x 4-matrix of unknowns. This ideal has n = 10, dim(/) = 7 
and homog(I) = 4. Hence T(I) is a two-dimensional polyhedral complex. 
We regard T(I) as the tropicalization of the secant variety of the Veronese 
threefold in F 9 , i.e., the variety of symmetric 4 x 4-matrices of rank < 2, 
Applying our Gf an implementation (see Example 15.4ft . we find that T(I) is 
a simplicial complex consisting of 75 triangles, 75 edges and 20 vertices. □ 

Our next problem concerns tropical bases. A finite set {fi, . . . , ft} is a 
tropical basis of I if (/i, . . . , f t ) = I and T(I) = T(/i) Pi • • • Pi T(f t ). 

Problem 2.6. Compute a tropical basis of a given ideal / C C[x]. 

A priori, it is not clear that every ideal / has a finite tropical basis, but 
we shall prove this below. First, here is one case where this is easy: 

Example 2.7. If / = (/) is a principal ideal, then {/} is a tropical basis. 

□ 

In it was claimed that any universal Grobner basis of I is a tropical 
basis. Unfortunately, this claim is false as the following example shows. 

Example 2.8. Let I be the intersection of the three linear ideals {x + y, z), 
{x + z,y), and (y + z,x) in C[x,y,z]. Then / contains the monomial xyz, 
so T(I) is empty. A minimal universal Grobner basis of / is 

U = {x + y + z, x 2 y + xy 2 , y 2 z + yz 2 , x 2 z + xz 2 }, 

and the intersection of the four corresponding tropical surfaces in R 3 is the 
line w\ = W2 = W3. Thus IA is not a tropical basis of /. □ 

We now prove that every ideal I C C[x] has a tropical basis. By Lemma r2.21 
one tropical basis of a non-homogeneous ideal / is the dehomogenization of 
a tropical basis for I. Hence we shall assume that / is a homogeneous ideal. 

Tropical bases can be constructed from the Grobner fan of / (see [Ej) |17j ) 
which is a complete finite rational polyhedral fan in IR n whose relatively open 
cones are in bijection with the distinct initial ideals of /. Two weight vectors 
w,w' € W 1 lie in the same relatively open cone of the Grobner fan of I if and 
only if in w (I) = in w i(I). The closure of this cell, denoted by C W (I), is called 
a Grobner cone of /. The n-dimensional Grobner cones are in bijection with 
the reduced Grobner bases, or equivalently, the monomial initial ideals of 
/. Every Grobner cone of / is a face of at least one n-dimensional Grobner 
cone of /. If in„,(7) is not a monomial ideal, then we can refine w to -< w by 
breaking ties in the partial order induced by w with a fixed term order -< 
on C[x]. Let Q^ W {I) denote the reduced Grobner basis of / with respect to 
~< w . The Grobner cone of G^ W (I), denoted by C^ W (I), is an n-dimensional 
Grobner cone that has C W (I) as a face. The tropical variety T(I) consists 
of all Grobner cones C W (I) such that in w (I) does not contain a monomial. 
From the description of T(I) as fl/e/ ^"(/) ^ is clear that T(I) is closed. 
Thus we deduce that T(I) is a closed subfan of the Grobner fan. This 
endows the tropical variety T(I) with the structure of a polyhedral fan. 
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Theorem 2.9. Every ideal I C C[x] has a tropical basis. 

Proof. Let T be any finite generating set of I which is not a tropical basis. 
Pick a Grobner cone C W (I) whose relative interior intersects C\f^jrT{f) non- 
trivially and whose initial ideal in w (I) contains a monomial x m . Compute 
the reduced Grobner basis Q^ W {I) for a refinement < w of w, and let h be 
the normal form of x m with respect to Q^ W (I). Let / := x m — h. Since the 
normal form of x m with respect to G^(m w (I)) = {in w (g) : g £ G-< W (I)} ^ s 

and h is the normal form of x m with respect to Q <w (I), every monomial 
occurring in h has higher w- weight than x m . Moreover, h depends only 
on the reduced Grobner basis G^ W (I) and is independent of the particular 
choice of w in C W (I). Hence for any w' in the relative interior of C W (I), 
we have x m = m w /(f). This implies that the polynomial / := x m — h is a 
witness for the cone C W (I) not being in the tropical variety T(I). 

We now add the witness / to the current basis T and repeat the process. 
Since the Grobner fan has only finitely many cones, this process will termi- 
nate after finitely many steps. It removes all cones of the Grobner fan which 
violate the condition for T to be a tropical basis. □ 

We next show that tropical bases can be very large even for linear ideals. 
Let / be the ideal in C[x] generated by d linear forms X]j=i a ij x j where 

1 = 1, . . . , d and ((%■) is an integer d x n matrix of rank d. The tropical 
variety T(I) depends only on the matroid associated with /, and it is known 
as the Bergman fan of that matroid. The results on the Bergman fan proved 
in p]E] imply that the circuits in / form a tropical basis. A circuit of / is a 
non-zero linear polynomial / E I of minimal support. The following result 
answers the question which was posed in |16| §5]. 

Theorem 2.10. For any 1 < d < n, there is a linear ideal I in C[xi, . . . , x n ] 
such that any tropical basis of linear forms in I has size at least n J d+l y) • 

Proof. Suppose that all d x d-minors of the coefficient matrix (a^) are non- 
zero. Equivalently, the matroid of / is uniform. There are ( n _^ +1 ) circuits 
in /, each supported on a different (n — d + l)-subset of {x\, . . . , x n }. Since 
the circuits form a tropical basis of / and each circuit has support of size 
n — d + 1, the tropical variety T(I) consists of all vectors w G M n whose 
smallest d + 1 components are equal. The latter condition is necessary and 
sufficient to ensure that no single variable in a circuit becomes the initial 
form of the circuit with respect to w. Consider any vector w € W 1 satisfying 

w h =w i2 = ■■■ = w id < minfwj : j 6 {1, . . . ,n}\{i 1 ,i 2 , ■ ■ -,id})- 

Since w any tropical basis of linear forms in / contains an / such that 

in w{f) S {xi x , . . . ,X{ d }. This implies that / is one of the d circuits whose 
support contains the n — d variables Xj with j {ii, . . . , id}. The support 
of each circuit has size n — d + 1, hence contains n — d + 1 distinct (n — d)- 
subsets. There are (^) (n — <i)-subsets of {xi, . . . , x n } to be covered. Hence 
any tropical basis consisting of linear forms has size at least n _ d+1 (^) • □ 
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Example 2.11. Let d = 3,n = 5. The Bergman fan T(J) corresponds to 
the line in tropical projective 4-space which consists of the five rays in the 
coordinate directions. We have n _^ +1 = 10/3. Hence this line is not a 
complete intersection of three tropical hyperplanes, but it requires four. □ 

3. Transversality and Connectivity 

In this section we assume that / is a prime ideal of dimension d in 
C[xi, . . . , x n ]. Then its tropical variety T(I) is called irreducible. It is a 
subfan of the Grobner fan of I and, by the Bieri-Groves Theorem |3|llS|. all 
facets of T(I) are cones of dimension d. A cone of dimension d— 1 in T(I) is 
called a ridge of the tropical variety T(I). A ridge path is a sequence of facets 
Fi,F2, . . . , Fk such that Fi n Fi + \ is a ridge for all i G {1, 2, . . . , k — 1}. Our 
objective is to prove the following result, which is crucial for the algorithms. 

Theorem 3.1. Any irreducible tropical variety T{I) is connected in codi- 
mension one, i.e., any two facets are connected by a ridge path. 

The proof of this theorem will be based on the following important lemma. 

Lemma 3.2. (Transverse Intersection Lemma) 

Let I and J be ideals in C[x\, . . . , x n ] whose tropical varieties T{I) and T{J) 
meet transversally at a point w 6 M. n . Then w £ T{I + J). 

By "meet transversely" we mean that if F and G are the cones of T(I) 
and 7~( J) which contain w in their relative interior, then MF + MG = M. n . 

This lemma implies that any transverse intersection of tropical varieties 
is a tropical variety. In particular, any transverse intersection of tropical 
hypersurfaces is a tropical variety, and such a tropical variety is defined by 
an ideal which is a complete intersection in the commutative algebra sense. 

Corollary 3.3. For any two ideals I and J in C[x\, . . . ,x n ] we have 

T(I + J) C T{I) n T(J). 

Equality holds if the latter intersection is transverse at every point except 
the origin and the two fans meet in at least one point other than the origin. 

Proof. We have T{I) n T(J) = f] feI T(f) n f] feJ T(f) = f\ feIyjJ T{f). 
Clearly, this contains T(I + J) = f] fEl+J T(f). If T(I) and T(J) intersect 
transversally and w is a point of T(I) PlT(J) other than the origin then the 
preceeding lemma tells us that w £ T{I + J). Thus T{I + J) contains every 
point of T{I) n T( J) except possibly the origin. In particular, T{I + J) 
is not empty. Every nonempty fan contains the origin, so we see that the 
origin is in T(I + J) as well. □ 

We first derive Theorem 13. II from Lemma 13.21 which will be proved later. 
We must at this point address an annoying technical detail. The subset 
T(I) C W 1 depends only on the ideal /Cb^ 1 ] generated by I in the Laurent 
polynomial ring C[x x , . . . , x^ 1 ]- (This is easy to see: if I\ and I2 generate the 



8 



T. BOGART, A. JENSEN, D. SPEYER, B. STURMFELS, R. THOMAS 



same ideal in C[x ] and w T{Ii) then there is a polynomial / 6 I\ such 
that in w (f) is a monomial. There is some monomial m such that mf £ -Z2, 
then m w (mf) is a monomial and to Tfa).) From a theoretical perspective 
then, it would be better to directly work with ideals in Cfx^ 1 ]. One reason 
is the availability the symmetry group GL n (Z) of the multiplicative group of 
monomials. The action of this group transforms T(I) by the obvious action 
on M n . This symmetry will prove invaluable for simplifying the arguments 
in this section. Therefore, in this section, we will work with ideals in Cbx:^ 1 ]. 
Computationally, however, it is much better to deal with ideals in C[x] as 
it is for such ideals that Grobner basis techniques have been developed and 
this is the approach we take in the rest of the paper. 

Note that, if I C C[x] is prime then so is the ideal it generates in Cfx 1 * 11 ]. 
We will signify an application of the GL n (Z) symmetry by the phrase "mak- 
ing a multiplicative change of variables" . The polyhedral structure on T{I) 
induced by the Grobner fan of I may change under a multiplicative change 
of variables of ICfx 1 * 1 ] in Cfx^ 1 ], but all of the properties of T(I) that are 
of interest to us depend only on the underlying point set. 

Proof of Theorem \3.1\ As discussed, we replace I by the ideal it generates 
in C[x ±:L ] and, by abuse of notation, continue to denote this ideal as I. The 
proof is by induction on d = dim(T(I)). If d < 1 then the statement is triv- 
ially true. We now explain why the result holds for d = 2. By a multiplica- 
tive change of coordinates, it suffices to check that T(I) n {x n = 1} is con- 
nected. Let K be the Puiseux series field over C. Let I' C K[x±, . . . , x n _i] 
be the prime ideal generated by I via the inclusion C[x n ] — > K. By Lemma 
11.11 the tropical variety of I' is T{I) n {x n = 1}. In 5 it was shown that 
the tropical variety of I' is connected whenever I 1 is prime. We conclude 
that T(I) n {x n = 1} is connected, so our result holds for d = 2. 

We now suppose that d > 3. Let F and F' be facets of T(I). We can find 

H = u n ) G R" : a\U\ -\ h a n u n = } 

such that ax, . . . , a n are relatively prime integers, both HnF and HCiF' are 
cones of dimension d—1, and H intersects every cone of T{I) except for the 
origin transversally. To see this, select rays w and w' in the relative interiors 
of F and F' . By perturbing w and w' slightly, we may arrange that the span 
of w and w' does not meet any ray of T(I) - here it is important that d > 3. 
Now, taking H to be the span of w, w' and a generic (n — 3)-plane, we get 
that H also does not contain any ray of T(I) and hence does not contain any 
positive dimensional face of T(I). So H is transverse to T(I) everywhere 
except at the origin. Since H n F and HnF' are positive-dimensional (as 
d > 2), the hyperplane H does intersect T(I) at points other than just the 
origin. The hyperplane H is the tropical hypersurface of a binomial, namely, 
H = T((/ u », where 

fu = n - n (w)~ a, > 

i:a;>0 j:aj<0 
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and u = (ui, U2, ■ ■ ■ ,u n ) is an arbitrary point in the algebraic torus (C*) n . 
Our transversality assumption regarding H and Lemma 13.21 imply that 

(1) h n T(i) = T((/„))nr(/) = r(j + (/„)). 

Since I is prime of dimension d, and f u G" I, the ideal I + {f u ) has dimension 
d— 1 by Krull's Principal Ideal Theorem j^l Theorem 10.1]. If 1+ (f u ) were 
a prime ideal then we would be done by induction. Indeed, this would 
imply that there is a ridge path between the facets H n F and H f] F' in 
the (d — l)-dimensional tropical variety (pQ). Since d > 3, the (c£ — 1)- and 
(d — 2)-dimensional faces of H Pi T(7) arise uniquely from the intersections 
of 77 with d- and (c£ — l)-dimensional faces of T{I). Hence this path is also 
a ridge path considered as a path in T(I). 

Let V(J) denote the subvariety of the algebraic torus (C*) n defined by 
an ideal J C C[xf 1 , . . . , x^ 1 ]. The tropical variety in (JTJ depends only on 
the subvariety of (C*) n defined by our ideal I + (f u )- This subvariety is 

(2) v(i+(f u )) = v(i)nv(f u ) = v(i) nu- 1 -v(h). 

Here 1 denotes the identity element of (C*) n . For generic choices of the 
group element u G (C*) n , the intersection © is an irreducible subvariety of 
dimension d — 1 in (C*) n . This follows from Kleiman's version of Bertini's 
Theorem |101 Theorem III. 10. 8], applied to the algebraic group (C*) n . Hence 
(^P) is indeed an irreducible tropical variety of dimension d — 1, defined by 
the prime ideal I + (f u )- This completes the proof by induction. □ 

Proof of Lemma \3.'A Again, we replace I C C[x] by the ideal it generates 
in Cfx^ 1 ]. Let F be the cone of T{t) which contains w in its relative 
interior and G the cone of T(J) which contains w in its relative interior. 
Our hypothesis is that F and G meet transversally at w, that is, 

RF + RG = R n . 

We claim that the ideal in„,(7) is homogeneous with respect to any weight 
vector v G RF or, equivalently (see Proposition 12. 4j) . that m v (in w (I)) = 
m w (I). According to Proposition 1.13 in ^7], for e a sufficiently small 
positive number, in TO 4. e „(J) = in v (m w (I)). The vector w + ev is in the 
relative interior of F so in w+ev (I) = in w (I). By the same argument, the 
ideal in w (J) is homogeneous with respect to any weight vector in RG. 

After a multiplicative change of variables in Cfx^ 1 , . . . , x^ 1 ] we may as- 
sume that w = ei, K{ei, e2, • • • , e s } C RF and K{ei, e s+ i, . . . , e n } C EG. 
We change the notation for the variables as follows: 

t = X!, y = (y 2 , ...,y s ) = (x 2 , . . -,x s ), z = (z s+1 , . . . , z n ) = (x s+ i, . . .,x n ). 

The homogeneity properties of the two initial ideals ensure that we can 
pick generators fi(z), ... , f a (z) for m w (I) and generators gi(y), . . -,gb(y) 
for in w (J). Since in„,(7) is not the unit ideal, the Laurent polynomials fi(z) 
have a common zero Z — (Z s -\-i 7 . . . , Z n ^j G (C*) n s , and likewise the Laurent 
polynomials gj(y) have a common zero Y = {Y 2 , ■■■ , Y s ) G (C*) s_1 . 
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Next we consider the following general chain of inclusions of ideals: 

(3) in w (I) -m w (J) C in™ (I • J) C in w (I n J) C in w (I) n m w (J). 

The product of two ideals which are generated by (Laurent) polynomials in 
disjoint sets of variables equals the intersection of the two ideals. Since the 
set of y- variables is disjoint from the set of z- variables, it follows that the 
first ideal in J2J) equals the last ideal in l|3*|). In particular, we conclude that 

(4) in™ (I n J) = in™ (I) n in™ (J). 
We next claim that 

(5) m w (I + J) = in™ (J) + in™ (J). 

The left hand side is an ideal which contains both in™ (I) and in™ (J), so it 
contains their sum. We must prove that the right hand side contains the left 
hand side. Consider any element / + g £ / + J where / £ I and g £ J. Let 
/ = fo(y,z) + t ■ fx(t,y,z) and g = g (y, z) + t ■ gi(t,y, z). We have the 
following representation for some integer o > and non-zero polynomial ho'- 

f + g = t a ■h (y,z)+t a+1 .h!(t,y,z). 

If a = then we conclude 

in™(/ + 5) = h {y,z) = f (y, z) + g (y, z) £ in™ (J) + in™( J). 

If a > 1 then Jq = —go lies in in™(/) n in™ (J). In view of (HJ), there exists 
p £ I n J with /o = -go = in™(p). Then / + g = (/ - p) + (g + p) and 
replacing / by (/ — p)/t and g by (g + p)/t puts us in the same situation 
as before, but with a reduced by 1. By induction on a, we conclude that 
in™(/ + g) is in in™(/) + in™ (J), and the claim © follows. 

For any constant T £ C*, the vector (T, Y%, . . . , Y s , Z s+ \, . . . , Z n ) is a 
common zero in (C*) n of the ideal (JSJ). We conclude that in™(/ + J) is not 
the unit ideal, so it contains no monomial, and hence w £ T{I + J). □ 

4. Algorithms 

In this section we describe algorithms for solving the computational prob- 
lems raised in Sectional The emphasis is on algorithms leading to a solution 
of Problem 12.31 for prime ideals, taking advantage of Theorem 13.11 Recall 
that we only need to consider the case of homogeneous ideals in C [x] . 

In order to state our algorithms we must first explain how polyhedral 
cones and polyhedral fans are represented. A polyhedral cone is represented 
by a canonical minimal set of inequalities and equations. Given arbitrary 
defining linear inequalities and equations, the task of bringing these to a 
canonical form involves linear programming. Representing a polyhedral fan 
requires a little thought. We are rarely interested in all faces of all cones. 

Definition 4.1. A set S of polyhedral cones in R n is said to represent a fan 
J- in lR n if the set of all faces of cones in S is exactly J-. 
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A representation may contain non-maximal cones, but each cone is repre- 
sented minimally by its canonical form. A Grobner cone C W {I) is represented 
by the pair (G^ w (m w (I)),Q^ w (I)) of marked reduced Grobner bases, where 
-< is some globally fixed term order. In a marked Grobner basis the initial 
terms are distinguished. The advantage of using marked Grobner bases is 
that the weight vector w need not be stored - we can deduce defining in- 
equalities for its cone from the marked reduced Grobner bases themselves, 
see Example 15.11 This is done as follows; see |17| proof of Proposition 2.3]: 

Lemma 4.2. Let I C C[x] be a homogeneous ideal, -< a term order and 
w G M. n a vector. For any other vector w' G M. n : 

w G C W (I) <S=^ V/ € Q-i w {I) ■ m w (m w <(f)) = in w (f). 

Our first two algorithms perform polyhedral computations, and they solve 
Problem By the support of a fan we mean the union of its cones. Recall 
that, for a polynomial /, the tropical hypersurface T(f) is the union of 
the normal cones of the edges of the Newton polytope New(/). The first 
algorithm computes these cones. 

Algorithm 4.3. Tropical Hypersurface 
Input: / G C[x]. 

Output: A representation S of a polyhedral fan whose support is T(f). 
{ 

5:=0; 

For every vertex v G New(/) 
{ 

Compute the normal cone C of v in New(/); 
S := Su{the facets of C}; 

} 

} 

Let T\ and be polyhedral fans in W 1 . Their common refinement is 

T\ A T 2 ■= {C\ n C2}(c 1 ,c 2 )eFixF2- 

To compute a common refinement we simply run through all pairs of cones 
in the fan representations and bring their intersection to canonical form. 
The canonical form makes it easy to remove duplicates. 

Algorithm 4.4. Common Refinement 

Input: Representations Si and 52 for polyhedral fans Ti and Ti- 
Output: A representation S for the common refinement Ti A Ti- 
{ 

S :=0; 

For every pair (Ci, C2) G Si x S2 
S :=Su{CinC 2 }; 

} 
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If refinements of more than two fans are needed, Algorithm 14.41 can be 
applied successively. Note that the intersection of the support of two fans 
is the support of the fans' common refinement. Hence Algorithm 14.41 can 
be used for computing intersections of tropical hyper surf aces. This solves 
Problem 12.11 but the output may be a highly redundant representation. 

Recall (from the proof of Theorem I2.9j) that a witness / G I is a poly- 
nomial which certifies T{f) n rel vat{C w {I)) = 0. Computing witnesses is 
essential for solving Problems 12.31 and 12.61 The first step of constructing 
a witness is to check if the ideal vn w (I) contains monomials, and, if so, 
compute one such monomial. The check for monomial containment can be 
implemented by saturating the ideal with respect to the product of the vari- 
ables (cf. |171 Lemma 12.1]). Knowing that the ideal contains a monomial, 
a simple way of finding one is to repeatedly reduce powers of the product of 
the variables by applying the division algorithm until the remainder is 0. 

Algorithm 4.5. Monomial in Ideal 

Input: A set of generators for an ideal / C C[x]. 

Output: A monomial m £ I if one exists, no otherwise. 

{ 

If ((/ : xi ■ ■ ■ / (1)) return no; 

TX\> • — X ^ ■ ■ ■ CCfi^ 

While (m £ I) m := m ■ x\ ■ ■ ■ x n ; 
Return m; 

} 

Remark 4.6. To pick the smallest monomial in / with respect to a term 
order, we first compute the largest monomial ideal contained in / using |14[ 
Algorithm 4.2.2] and then pick the smallest monomial generator of this ideal. 

Constructing a witness from a monomial was already explained in the 
proof of Theorem 12. 91 We only state the input and output of this algorithm. 

Algorithm 4.7. Witness 

Input: A set of generators for an ideal / C C[x] and a vector w G M n with 
in w (I) containing a monomial. 

Output: A polynomial / 6 I such that the tropical hypersurface T(f) and 
the relative interior of C W (I) have empty intersection. 

Combining Algorithm l4 . 51 and Algorithm l4. 7l with known methods (e.g. [171 
Algorithm 3.6]) for computing Grobner fans, we can now compute the trop- 
ical variety T(I) and a tropical basis of /. This solves Problem 12.31 and 
Problem 12.61 This approach is not at all practical, as shown in Section 6. 

We will present a practical algorithm for computing T{I) when / is prime. 
An ideal / C C[x] is said to define a tropical curve if dim(/) = l + homog(/). 
Our problems are easier in this case because a tropical curve consists of only 
finitely many rays and the origin modulo the homogeneity space. 

Algorithm 4.8. Tropical Basis of a Curve 

Input: A set of generators Q for an ideal / defining a tropical curve. 
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Output: A tropical basis Q' of I. 
{ 

Compute a representation S of f\ g& g T(g); 
For every C £ S 

{ 

Let w be a generic relative interior point in C; 
If (in w (I) contains a monomial) 

then add a witness to Q and restart the algorithm; 

} 

Q' := G; 

} 

Proof of correctness. The algorithm terminates because I has only finitely 
many initial ideals and at least one is excluded in every iteration. If a vector 
w passes the monomial test (which verifies w € T(I)) then C has dimension 
or 1 modulo the homogeneity space since we are looking at a curve and 
w is generic in C. Any other relative interior point of C would also have 
passed the monomial test. (This property fails in higher dimensions, when 
T(I) is no longer a tropical curve). Hence, when we terminate only points 
in the tropical variety are covered by S. Thus Q' is a tropical basis. □ 

In the curve case, combining Algorithms 14.31 and 14.41 with Algorithm 14.81 
we get a reasonable method for solving Problem [231 This method is used as 
a subroutine in Algorithm 14.101 below. In the remainder of this section we 
concentrate on providing a better algorithm for Problem 12.31 in the case of a 
prime ideal. The idea is to use connectivity to traverse the tropical variety. 

The next algorithm is an important subroutine for us. We only specify 
the input and output. This algorithm is one step in the Grobner walk 0j. 

Algorithm 4.9. Lift 

Input: Marked reduced Grobner bases Q^{I) and G^ w (in w (I)) where w £ 
C-ji(I) is an unspecified vector and -< and -<' are unspecified term orders. 
Output: The marked reduced Grobner basis Q <w {I). 

We now suppose that / is a monomial-free prime ideal with d = dim(J), 
and -< is a globally fixed term order. We first describe the local computations 
needed for a traversal of the d-dimensional Grobner cones contained in T(I). 

Algorithm 4.10. Neighbors 

Input: A pair (G^ w (m w (I)),Q^ w (^)) such that in w (I) is monomial-free and 
C W {I) has dimension d. 

Output: The collection N of pairs of the form ,(m w '(I)),G^ ,{!)) 
where one vJ is taken from the relative interior of each (i-dimensional Grobner 
cone contained in T{I) that has a facet in common with C W (I). 
{ 

iV:=0; 

Compute the set T of facets of C W (I); 
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Figure 1. A projective drawing of the situation in Algo- 
rithm |4~TU1 with T{I) on the left and T(in u (/)) on the right. 

For each facet F £ J- 
{ 

Compute the initial ideal J := in u (J) 

where u is a relative interior point in F; 
Use Algorithm 14,81 and Algorithm 14,41 to produce a relative 

interior point v of each ray in the curve T(J); 
For each such v 

{ 

Compute (<^ v (in v (J)),^ v (J)) = (^ Vu (in v ( J)), £^ Vu ( J)); 
Apply Algorithm EH1 to Q <W {I) and ^ Vu (J) to get G^JI); 
iV:=iVU{(^ Vu (in v (J)),^ Vu (/))}; 

} 

} 

} 

Proof of correctness. Facets and relative interior points are computed 
using linear programming. Figure ^ illustrates the choices of vectors in 
the algorithm. The initial ideal in u (7) is homogeneous with respect to the 
span of F. Hence its homogeneity space has dimension d — 1. The Krull 
dimension of C[x]/in u (7) is d. Hence in u (I) defines a curve and T(in u (/)) 
can be computed using Algorithm 14.81 The identity in v (in u (/)) = in u+EV (7) 
for small e > 0, see |17| Proposition 1.13], implies that we run through all 
the desired m. w i (I) where vJ = u + ev for small e > 0. The lifting step can 
be carried out since u G C^ W (I). □ 

Algorithm 4.11. Traversal of an Irreducible Tropical Variety 

Input: A pair (Q^ w (m w (I)),Q^ w (I)) such that in to (J) is monomial free and 

C W (I) has dimension d. 

Output: The collection T of pairs of the form (Q^ , (in^/(J)), , (I)) 
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where one w' is taken from the relative interior of each d-dimensional Grobner 

cone contained in T(I). The union of all the C w i(I) is T(I). 

{ 

T-={{Q^mv{I)),g <w (I))}\ 
Old : = 0; 

While (T ^ Old) 
{ 

Old := T; 

T := TU Neighbors(T); 

} 

} 

Proof of correctness. By Neighbors (T) we mean the union of all the output 
of Algorithm I4.1UI applied to all pairs in T. The algorithm computes the 
connected component of the starting pair. Since I is a prime ideal, Theorem 
13. II implies that the union of all the computed C w i(I) is T{I). □ 
To use Algorithm 14.111 we must know a starting <i-dimensional Grobner 
cone contained in the tropical variety. One inefficient method for finding 
one would be to compute the entire Grobner fan. Instead we currently use 
heuristics, which are based on the following probabilistic recursive algorithm: 

Algorithm 4.12. Starting Cone 

Input: A marked reduced Grobner basis Q for an ideal I whose tropical 
variety is pure of dimension d = dim(/). A term order -< for tie-breaking. 
Output: Two marked reduced Grobner bases: 

• One for an initial ideal m w /(I) without monomials, where the homo- 
geneity space of m w /(I) has dimension d. The term order is -< w i. 

• A marked reduced Grobner basis for / with respect to -< w '. 

{ 

If (dim(7) = homog(/)) 

Return (^(J), £<(/)); 
If not 
{ 

Repeat 
{ 

Compute a random reduced Grobner basis of /; 
Compute a random extreme ray w of its Grobner cone; 

} 

Until (in„,(J) is monomial free); 
Compute g~> w (I); 

{Qlnit,QFull)-= Starting Cone(£^ ffi (in w (I))); 
Apply Algorithm to Q^ W {I) and SfuII 

to get a marked reduced Grobner basis Q' for /; 
Return {Qi n it,Q')\ 

} 

} 
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5. Software and Examples 

We implemented the algorithms of Section |1] in the software package 
Gfan 9 . Gfan uses the library cddlib [Jj for polyhedral computations such 
as finding facets and extreme rays of cones and bringing cones to canonical 
form. We require our ideals to be generated by polynomials in Q[x]. Exact 
arithmetic is done with the library gmp (Sj. This is needed both for poly- 
hedral computations and for efficient arithmetic in Q[x]. In this section we 
illustrate the use of Gfan in computing various tropical varieties. 

Example 5.1. We consider the prime ideal / C C[a, b, c, d, e, f, g] which is 
generated by the 3x3 minors of the generic Hankel matrix of size 4x4: 
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Its tropical variety is a 4-dimensional fan in M7 with 2-dimensional homo- 
geneity space. Its combinatorics is given by the graph in Figure [3 To 
compute T(I) in Gfan, we write the ideal generators on a file hankel . in: 

°/ more hankel. in 

{-c"3+2*b*c*d-a*d~2-b~2*e+a*c*e ,-c~2*d+b*d~2+b*c*e-a*d*e-b~2*f +a*c*f , 
-c*d~2+c~2*e+b*d*e-a*e~2-b*c*f +a*d*f , -d"3+2*c*d*e-b*e~2-c"2*f +b*d*f , 
-c~2*d+b*d~2+b*c*e-a*d*e-b~2*f +a*c*f , -c*d~2+2*b*d*e-a*e~2-b"2*g+a*c*g, 
-d~3+c*d*e+b*d*f -a*e*f -b*c*g+a*d*g, -d~2*e+c*e~2+c*d*f -b*e*f -c"2*g+b*d*g, 
-c*d~2+c~2*e+b*d*e-a*e~2-b*c*f +a*d*f , -d~3+c*d*e+b*d*f -a*e*f -b*c*g+a*d*g, 
-d~2*e+2*c*d*f-a*f ~2-c~2*g+a*e*g,-d*e"2+d~2*f +c*e*f -b*f ~2-c*d*g+b*e*g, 
-d~3+2*c*d*e-b*e~2-c~2*f +b*d*f , -d"2*e+c*e~2+c*d*f-b*e*f -c"2*g+b*d*g, 
-d*e~2+d~2*f +c*e*f-b*f ~2-c*d*g+b*e*g,-e~3+2*d*e*f-c*f ~2-d"2*g+c*e*g> 

We then run the command 
gf an_tropicalstartingcone < hankel. in > hankel. start 

which applies Algorithm 14.121 to produce a pair of marked Grobner bases. 
This represents a maximal cone in T(J), as explained prior to Lemma 14.21 

"L more hankel . start 
{ 

c*f "2-c*e*g , 
b*f "2-b*e*g, 
b*e*f +c"2*g, 
b*e*2+c"2*f , 
b"2*g-a*c*g, 
b"2*f-a*c*f , 
b"2*e-a*c*e , 
a*f "2-a*e*g, 
a*e*f +b*c*g , 
a*e~2+b*c*f } 
{ 

c*f "2+e"3-2d*e*f +d~2*g-c*e*g, 
b*f "2+d*e"2-d"2*f -c*e*f +c*d*g-b*e*g, 
b*e*f +d"2*e-c*e"2-c*d*f +c"2*g-b*d*g, 
b*e"2+d"3-2c*d*e+c~2*f -b*d*f , 
b"2*g+c"2*e-b*d*e-b*c*f +a*d*f -a*c*g, 
b"2*f +c"2*d-b*d"2-b*c*e+a*d*e-a*c*f , 
b"2*e+c"3-2b*c*d+a*d~2-a*c*e ) 
a*f "2+d"2*e-2c*d*f +c~2*g-a*e*g, 
a*e*f +d"3-c*d*e-b*d*f +b*c*g-a*d*g , 
a*e"2+c*d"2-c"2*e-b*d*e+b*c*f -a*d*f } 



COMPUTING TROPICAL VARIETIES 



17 



Using Lemma 14.21 we can easily read off the canonical equations and equal- 
ities for the corresponding Grdbner cone C W {I). For example, the polyno- 
mials c/ 2 — ceg and c/ 2 + e 3 — 2def + d 2 g — ceg represent the equation 

w c + 2wf = w c + w e + Wg 

and the inequalities 

w c + 2wf < min{3w e , Wd + w e + Wf,2w,i + w g ,w c + w e + w g } . 
At this point, we could run Algorithm 14. 1 II using the following command: 
gf an_tropicaltraverse < hankel. start > hankel.out 

However, we can save computing time and get a better idea of the struc- 
ture of T(I) by instructing Gf an to take advantage of symmetries of / as it 
produces cones. The only symmetries that can be used in Gfan are those 
that simply permute variables. The output will show which cones of T(I) 
lie in the same orbit under the action of the symmetry group we provide. 

Our ideal / is invariant under reflecting the 4 x 4-matrix along the anti- 
diagonal. This reverses the variables a,b,...,g. To specify this permutation, 
we add the following line to the bottom of the file hankel . start: 

{(6,5,4,3,2,1,0)} 

We can add more symmetries by listing them one after another, separated 
by commas, inside the curly braces. Gfan will compute and use the group 
generated by the set of permutations we provide, and it will return an error 
if we input any permutation which does not keep the ideal invariant. 

After adding the symmetries, we run the command 
gf an_tropicaltraverse — symmetry < hankel. start > hankel.out 
to compute the tropical variety. We show the output with some annotations: 



X more hankel.out 



A short list of basic data: the dimensions 
of the ambient space, of T(I), and of 
its homogeneity space, and also the face 
numbers (/-vector) of T(I) and the order 
of symmetry group specified in the input. 
A basis for the homogeneity space. The 
rays are considered in the quotient of M7 
modulo this 2-dimensional subspace. 



Ambient dimension: 7 



Dimension of homogeneity space: 2 
Dimension of tropical variety: 4 
Simplicial: true 



Order of input symmetry group: 2 
F-vector: (16,28) 



Modulo the homogeneity space: 
{(6,5,4,3,2,-1,0) , 
(6,4,3,2,1,0,-1)} 



Rays : 
{0: ( 



(-1,0,0,0,0,0,0) , 
(-5,-4,-3,-2,-1,0,0) , 



2: 
3: 
4: 
5: 
6: 
7: 
8: 
9: 
10 
11 
12 
13 
14 
15 



(1,0,0,0,0,0,0) , 
(5,4,3,2,1,0,0), 
(2,1,0,0,0,0,0) , 
(4,3,2,1,0,0,0), 
(0,-1,0,0,0,0,0) , 
(6,5,4,3,2,0,0) , 
(3,2,1,0,0,0,0), 
(0,0,-1,0,0,0,0) , 
(0,0,0,0,-1,0,0) 
(0,0,0,-1,0,0,0) 



(-6,-4,-3,-3,-1,0,0) , 
(-3,-2,-2,-1,-1,0,0) , 



The direction vectors of the tropical 
rays. Since the homogeneity space is 
positive-dimensional, the directions are 
not uniquely specified. For instance, 
the vectors (—5, —4, —3, —2, —1, 0, 0) and 
(0, 0, 0, 0, 0, 0, —1) represent the same ray. 
Note that Gfan uses negated weight vec- 
tors. 



(3,2,2,1,1,0,0) , 
(3,2,2,0,1,0,0)} 
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Figure 2 . The tropical variety of the ideal generated by the 
3x3 minors of the generic 4x4 Hankel matrix. 

Rays incident to each 
dimension 2 cone: 
{{2,6}, {3,7}, 
{2,4}, {3,5}, 
{4,9}, {5,10}, 
{4,8}, {5,8}, 
{8,11}, 

{0,12}, {1,12}, 
{0,1}, 

{1,6}, {0,7}, 
{1,9}, {0,10}, 
{0,13}, {1,13}, 
{6,14}, {7,14}, 
{9,13}, {10,13} 
{6,10}, {7,9}, 
{6,7}, 
{11,12}, 
{11,15}, 
{14,15}} 

The further output, which is not displayed here, shows that the 16 rays 
break down into 5 orbits of size 2 and 6 orbits of size 1. 

Using the same procedure, we now compute several more examples. 

Example 5.2. Let / be the ideal generated by the 3x3 minors of the 
generic 5x5 Hankel matrix. We again use the symmetry group Z/2. The 
tropical variety is a graph with vertex degrees ranging from 2 to 7. 

Ambient dimension: 9 

Dimension of homogeneity space: 2 

Dimension of tropical variety: 4 

Simplicial: true 

F-vector: (28,53) 

Example 5.3. Let I be the ideal generated by the 3x3 minors of a generic 
3x5 matrix. We use the symmetry group S§ x S3, where £5 acts by per- 
muting the columns and 53 by permuting the rows. 

Ambient dimension: 15 

Dimension of homogeneity space: 7 

Dimension of tropical variety: 12 

Simplicial: true 

F-vector: (45,315,930,1260,630) 



The cones in T(I) are listed from highest 
to lowest dimension. Each cone is named 
by the set of rays on it. There are 28 two- 
dimensional cones, broken down into 11 
orbits of size 2 and 6 orbits of size 1. 
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Example 5.4. Let I be the ideal generated by the 3x3 minors of a generic 
4x4 symmetric matrix. We use the symmetry group £4 which acts by 
simultaneously permuting the rows and the columns. 

Ambient dimension: 10 
Dimension of homogeneity space: 4 
Dimension of tropical variety: 7 
Simplicial: true 
F-vector: (20,75,75) 

If we take the 3x3 minors of a generic 5x5 symmetric matrix then we get 

Ambient dimension: 15 
Dimension of homogeneity space: 5 
Dimension of tropical variety: 9 
Simplicial: true 

F-vector: (75, 495, 1155, 855) 

Example 5.5. Let I be the prime ideal of a pair of commuting 2x2 matrices. 
That is, I C C[a, b, . . . , h] is defined by the matrix equation 

00 0-0 00 0= a 

The tropical variety is the graph K4, which Gf an reports as follows: 

Ambient dimension: 8 

Dimension of homeogeneity space : 4 

Dimension of tropical variety: 6 

Simplicial: true 

F-vector: (4,6) 

If / is the ideal of 3 x 3 commuting symmetric matrices then we get: 

Ambient dimension: 12 
Dimension of homeogeneity space : 2 
Dimension of tropical variety: 9 
Simplicial: false 

F-vector : (66 , 705 , 3246 , 7932 , 10888 , 8184 , 2745) 

6. Tropical variety versus Grobner fan 

In this paper we developed tools for computing the tropical variety T(I) 
of a d-dimensional homogeneous prime ideal I in a polynomial ring C[x]. We 
took advantage of the fact that, since / is homogeneous, the set T{I) has 
naturally the structure of a polyhedral fan, namely, T(I) is the collection 
of all cones in the Grobner fan of / whose corresponding initial ideal is 
monomial-free. A naive algorithm would be to compute the Grobner fan of 
/ and then retain only those d-dimensional cones who survive the monomial 
test (Algorithm 14.5(1 . The software Gfan also computes the full Grobner fan 
of /, and so we tested this naive algorithm. We found it to be too inefficient. 
The reason is that the vast majority of d-dimensional cones in the Grobner 
fan of I are typically not in the tropical variety T(I). 

Example 6.1. Consider the ideal / in Example 15.11 which is generated by 
the 3 x 3-minors of a generic 4 x 4-Hankel matrix. Let J = in w (I) be its 
initial ideal with respect to the first vector w in the list of rays. The initial 
ideal J defines a tropical curve consisting of five rays and the origin. The 
curve is a subfan of the much more complicated Grobner fan of J. The 
Grobner fan is full-dimensional in W. 7 with Cq(J) being three-dimensional. 
Its f-vector equals (1, 7167, 32656, 45072, 19583). Of the 7167 rays only 5 are 
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in the tropical variety. The Grobner fan of J is the link of the Grobner fan 
of I at w. We were unable to compute the full Grobner fan of /. 

Example 6.2. Toric Ideals. Let / = (x u — x v : ^4u = ^4v) be the toric 
ideal of a matrix A £ 2^ xn Q f rank d. The ideal I is a prime of dimension d. 
The tropical variety T{I) coincides with the homogeneity space Cq{I) which 
is just the row space of A. Hence T(I) modulo C${I) is a single point. Yet, 
the Grobner fan of I can be very complicated, as it encodes the sensitivity 
information for an infinite family of integer programs |17| Chapter 7]. 

We next exhibit a family of ideals such that the number of rays in T{T) 
is constant while the number of rays in the Grobner fan of I grows linearly. 

Theorem 6.3. Fix n = 3,d=l and for any integer p > 1 consider the ideal 

I P = (x-(z + ir+ 2 ,y-(z-iy). 

Then T(I p ) consists of 4 rays but the Grobner fan of I p has > \{p+^) rays. 

Sketch of proof: The ideal I p is prime. Its variety is the parametric curve 
z i ► ((z+l) p+2 , [z— z) . The poles and zeros of this map are 0, — 1, +1, oo. 
The tropical variety of I p consists of the four rays defined by the valuations 
at these points. These rays are generated by the columns of 

p + 2 -p - 2 \ 
Op — p 
10 -1 J 

We examine the Grobner fan around the ray w = —(p + 2,p, 1). The initial 
ideal in^(I p ) equals the toric ideal (x— z p+2 , y—z p ). To see this, we note that 
the two generators of I p form a Grobner basis with respect to the underlined 
leading terms and in w (I p ) is generated by in w (g) for each g in this Grobner 
basis since w lies in this Grobner cone. The Grobner fan of m w (I p ) is the 
link at w of the Grobner fan of I p . To prove the theorem we show that the 
Grobner fan of in W) (/ p ) has at least |(p + 1) distinct Grobner cones. This 
implies, by Euler's formula, that the Grobner fan of in^(I p ) has at least 
\{p+ 1) rays and hence so does the Grobner fan of I p . 

To argue that the Grobner fan of m w (I p ) has at least |(p + 1) distinct 
Grobner cones we use the methods in |17j . More specifically, this involves 
first showing that the binomials gj := y J — z p ~ 2 ^~ "^a^ -1 for j = 1, . . . , Pj ^- 
are all in the universal Grobner basis of in w (I p ). Each monomial in a bino- 
mial in the universal Grobner basis of a toric ideal contributes a minimal 
generator to some initial ideal of the toric ideal. Thus there exist reduced 
Grobner bases of in w (I p ) in which the binomials gj are elements with leading 
term if for j = 1, . . . , £t-. This implies that these reduced Grobner bases 
are all distinct, which completes the proof. 

□ 

While the Grobner fan is a fundamental object which has had a range 
of applications (the Grobner walk 0j, integer programming (Example E2J)), 
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many computer algebra experts do not like it. Their view is that the Grobner 
fan is a combinatorial artifact which is marginal to the real goal of computing 
the variety of I. While this opinion has some merit, the story is entirely 
different for the subfan T{I) of the Grobner fan. In our view, the tropical 
variety is the variety of I. Every point on T{I) furnishes the starting system 
for a numerical homotopy towards the complex variety of I, see |18( Chapter 
3]. Thus computing T(I) is not only much more efficient than computing 
the Grobner fan of /, it is also geometrically more meaningful. 
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